A SUPERFLUID GYROSCOPE WITH COLD ATOMIC GASES 



S. Stringari 

Dipartimento di Fisica, Universita di Trento, 
and Istituto Nazionale per la Fisica della Materia, 
1-38050 Povo, Italy 
(February 1, 2008) 



Abstract 

A trapped Bose-Einstein condensed atomic gas containing a quantized vortex 
is predicted to exhibit precession after a sudden rotation of the confining 
potential. The equations describing the motion of the condensate are derived 
and the effects of superfluidity exphcitly pointed out. The dependence of the 
precession frequency on the relevant parameters of the problem is discussed. 
The proposed gyroscope is well suited to explore rotational effects at the level 
of single quanta of circulation. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj, 67.40.-w 



Typeset using REVTgX 



After the realization of Bose-Einstein condensation [|l| in magnetically trapped vapours 
of alkali atoms, cooled down to extremely low temperatures, the superfluid behaviour of 
these systems has become the object of extensive experimental work. This includes the 
study of rotational properties, like quantized vortices [^|[ and the quenching of the moment 
of inertia [Q, as well as the reduction of dissipative effects 0. 

The purpose of this work is to show that Bose-Einstein condensed (BEC) gases can be 
used to realize a quantum gyroscope where the effects of superfluidity show up in a very 
peculiar way. Superfluid gyroscopes have been already realized with liquid ^He and ^He 
and have been mainly used to investigate the nature of persistent currents in toroidal 
geometries. Most of experiments with helium gyroscopes operate with many quanta of 
circulation. Compared to liquid helium, trapped BEC gases are mesoscopic systems in the 
sense that the healing length, which provides a typical range of dynamic correlations and 
fixes the size of the vortex core, is smaller, but not extremely smaller than the size of the 
sample ^j. For the same reason the angular momentum Nh carried by a single quantized 
vortex can have visible effects on the global motion of the condensate as we will prove in 
this letter. A further feature that characterizes these systems is the very peculiar type of 
confinement which yields new possibilities for exploring superfluid phenomena. Important 
gyroscopic effects associated with the occurrence of vortex lines have been already observed 
in trapped BEC gases ||9|-[TT|. The authors of have succeeded in testing the quantization of 



the angular momentum of a single vortex line. To this purpose one generates a quadrupole 
deformation in the plane orthogonal to the vortex axis. The observed precession of the 
deformation is proportional to the angular momentum carried by the vortex, in accordance 
with the predictions of theory |[T^JT3|1 . Notice that in this experiment the vortex line is not 



affected by the precession. The authors of |10| and [11 1 have instead observed the precession 



of a vortex line either displaced or tilted from the symmetry axis of the condensate. In this 
case the motion of the condensate is not affected at a macroscopic level and the precession 
involves the change of the density on a more microscopic scale, of the order of the size of 
the vortex core. 



In the present work we discuss the macroscopic precession of the symmetry axis of a 
deformed condensate caused by the sudden rotation of the trap, in the presence of a quantized 
vortex hne. This precession corresponds to a full rotation of the condensate in 3D, which 
preserves the intrinsic shape of the system. 

The geometry of the proposed gyroscope is illustrated in fig.l. The system consists of a 
dilute and cold gas of atoms with mass m confined by an axi-symmetric trap of harmonic 
shape: 



where uj±_ and ujz are, respectively, the radial and axial trapping frequencies. The trapping 
potential is characterized by the deformation parameter 



which will be taken different from zero. As a consequence of this anisotropy the atomic 
cloud is also deformed. If the atoms are bosons, at very low temperatures they undergo 
a phase transition associated with the macroscopic occupation of a single-particle state 
(Bose- Einstein condensation). This happens for temperatures below the critical value kTc = 
0.94:huJhoN^^^ where Uho = {^z^±y^^ geometrical average of the trapping frequencies 

and N is the number of atoms. Typical values of Tc range between 0.1 and IfiK. Below Tc 
the system exhibits unique features which have been the object of systematic experimental 
and theoretical investigation. For our purposes the system should be very cold in order to 
minimize the dissipative effects produced by the thermal cloud which tend to destabilize 
the vortex line p!5|JT6|] . For this reason we will discuss here the behaviour of the gas at zero 
temperature. 

Let us consider a Bose-Einstein condensate in equilibrium in the trap, with no angular 
momentum, and let us suppose that the symmetry axis of the trapping potential is suddenly 
rotated in the xz plane through an angle ^o- After the rotation the gas is no longer in 
equilibrium and starts oscillating around the new symmetry axis. If the angle of rotation 
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is small the gas will keep its intrinsic shape and the motion will correspond to a periodic 
oscillation characterized by the angle 0{t), giving rise to the so-called scissors mode [p!^ . 



The motion is very different from the one of a classical body because the inertia is strongly 
suppressed by superfluidity. In particular the oscillation frequency does not vanish when 
the deformation (^) of the trap becomes small, as one would expect for a classical system. 

The calculation of the frequency of the scissors mode requires the solution of a non trivial 
many-body problem where the interactions among particles and the effects of quantum 
statistics play an important role. An exact solution of the problem can be obtained for 
large samples, such that the dimensionless parameter Na/aho is much larger than unity. 
This limit corresponds to the so-called Thomas-Fermi (TF) regime and is well achieved 
in many available experimental configurations. Here a is the s-wave scattering length (of 



the order of 10^^ — 10^^ fim) while aho = \Jf>'l ijn^ho) is the oscillator length fixed by the 
geometrical average of the trapping frequencies (of the order of a few microns in standard 
magnetic traps). In the TF limit, the equilibrium density takes the form of an inverted 
parabola: n(r) = m[/i — V"e2;t(r)]/(47r/i^a) for /i > Vext and zero elsewhere, where jj, = 
(c(j/jo/2)(15A^a/a/jo)^/^ is the chemical potential fixed by the normalization of the density. 
The surface where the density vanishes defines an ellipsoid whose principal radii R± and 
Rz are determined by the relationship uo']_R\ = uj'^Rl = 2/i/m. In the Thomas- Fermi limit 
the equations of motion take the simplified form of the hydrodynamic theory of superfiuids 
r?| . With respect to the equations of classical hydrodynamics they are characterized but 



the additional constraint of irrotationality. It is worth noticing that the "hydrodynamic" 
form of these equations is not the result of collisional processes as happens in classical gases, 
but is the consequence of superfluidity. These equations admit an analytic solution for the 
scissors mode with frequency 



Us 



This mode has recently been observed [Q confirming the predictions of theory with high 
accuracy. 



Let us discuss now the case in which the condensate has initially an intrinsic angular 
momentum. This can be due to a quantized vortex aligned along the symmetry axis of the 
condensate, carrying one unit h of angular momentum per particle. The size of the vortex 
core is fixed by the healing length ^ = (2m/i//;,^)~^/^ and is always smaller than the size 
of the condensate. Also in the presence of the vortex the sudden rotation of the trap will 
excite the scissors mode, but now the oscillation will precess due to the torque 

^ < L >= < [L, i/] >= -m < r X WVext > (4) 

produced by the external anisotropic potential (|I|). Here L is the angular momentum op- 
erator and the expectation value is taken on the quantum mechanical state of the system. 
The value of the torque along the y axis is d < Ly > /dt = —m{ijj_ — uj1)N < xz > and 
differs from zero immediately after the sudden rotation of the trap in the a;2;-plane, thereby 
causing the precession. 

In order to derive the equation for the rotation of the condensate we notice that, if the 
angle of rotation 60 is small compared to the deformation (0) of the trap, the shape of the 
sample is preserved. Under this condition the whole motion can be described in terms of 
two angles: the inclination angle 6 and the azimuthal angle (p. This is a remarkable feature 
because gases are highly compressible and their shape can change very easily. The values of 

6 and are related to the averages < xz > and < yz > hj the geometrical relations 

z?2 , ,2 , ,2 p2 , ,2 , ,2 
< XZ >= — — ^6'cos0, < yz >= — — -6sm(j) (5) 

7 UJz^± ' ^z^l. 



where R = \jR±Rz and we have assumed 6* -C e. 

The crucial problem now is to obtain the equations for the averages < xz > and < yz > 
in the presence of the vortex. This problem was solved in [|12| where it was shown that, in the 
Thomas- Fermi limit, the quadrupole operators q± = Y.k=ii^k ^Wk)zk excite, respectively, 
just one collective mode | ± > with frequency u!±. The solution of the Schrodinger equation, 
in the presence of a small perturbation generated by these operators, then takes the form 

I ^(t) >= e-^^o*/^ [1 > +c+e-^"+* I + > +c_e-^"-* | - >! (6) 



where | > is the equihbrium configuration, | ± > are eigenstates of the Hamiltonian 
corresponding to the elementary excitations created by q± and the complex parameters c± 
characterize the initial state of the system. One can also show that the strengths of the 
quadrupole operators q± are equal: < | q^q+ \ >=< | | >, while the splitting 

between the two frequencies takes the simple form Wl\ 



<lz> 

uj+-uj- = — — - — (7) 

m < x"^ + > 

In eq.(0) the quantity < > is the angular momentum per particle in the equilibrium 
configuration, while the average square radii can be easily evaluated in the Thomas-Fermi 
limit where one finds < >=< y"^ >= {2/7)fi/mu\ and < z"^ >= {2/7)^/ mul. Result 
(0) holds if the splitting is small compared to the unperturbed value @ of the collective 
frequency. 

The scenario emerging from the above discussion is quite clear. In the absence of vortices 
(< Iz >= 0) the two modes | ± > are degenerate and their frequency is given by (0). This is 
a simple consequence of time reversal symmetry. In the presence of the vortex there is a lift 
of degeneracy which is determined by the angular momentum of the system according to (|^ 
|r8|. By expressing the average values < xz > and < yz > in terms of the expectation values 



of the operators q± and using explicitly eq. ®, one derives the following time dependence 
for the angles 6 and 0: 

9(t) cos (pit) = a cos cj+t + a' sin u+t + P cos U-t + P' sin U-t (8) 



9(t) sin = —a' cos u+t + a sin u+t + [3' cos u^t — {3 sin uj^t (9) 

where the real coefficients a, a', (3 and (3' depend on the initial conditions for 9 and (j). For 
example the normal modes | ± > can be separately excited by choosing a = 9q, a' = (3 = 
f3' = and (3 = 9o, a = a' = (3' = respectively. In both cases the inclination angle 9 
remains constant, 9{t) = 9q, while the azimuthal angle precesses according to the laws u^t 
and —U-t, respectively. 
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In the most relevant case of a sudden rotation in the xz plane, the initial conditions 
instead correspond to 6{0) = 9q, 0'{O) = 0(0) = <f)'{0) = 0. In this case the parameters of 
eqs.(|[]^) are a = (3 = 9o/2, a' = j3' = and the solutions take the form 

Mt) = ii^i^ (10) 

and 

{uj+ + UJ^)t 



e{t) = ^ocos 



(11) 



The simultaneous excitation of the | ± > modes is easily understood by noticing that 
the rotation of the trap produces a change in the external potential of the form 6Vext = 
m{uj'j_ — uj1)9qxz. This gives rise to a perturbative term in the Hamiltonian, proportional to 
the combination {q+ + qJ) of the quadrupole operators introduced above. Equations ([T0|) and 
(p!T|) show that the scissors oscillation, characterized by the frequency {uj^+ujJ) /2 ~ ujsdssors, 
undergoes a precession dcp/dt fixed by the splitting (^. The relative precession can be 
expressed in the form 

7<4> f^^^^y. ^^^^ 



CU+ + cu„ 2 n (1 + A2)3/2 V aho 



where A = u!z/uj± = y (1 — e)/(l + e) and we have used the Thomas-Fermi result for < > 
and < >. The precession frequency d(f)/dt, and hence the ratio (0), depends explictly 
on the value of the angular momentum per particle which, in the case of a single quantized 
vortex aligned along the symmetry axis of the condensate, is given by < >= h. The ratio 
(|l^) depends also on the shape of the trap. This provides further flexibility to optimize 
the visibility of the precession. In typical experimental configurations the Thomas-Fermi 
combination (15A^a/a/jo)^^/^ is of the order of 10~^ so that for a highly elongated trap 



(A ^ 1) the relative precession is small. For values of A closer to unity the ratio (|T2|) 
becomes larger. With suitable choices of the parameters of the trap the precession can be 
easily of the order of 1 Hz and should be consequently observable by imaging the atomic 
cloud at different times. 



In conclusion we have shown that Bose-Einstein condensed gases confined in harmonic 
traps can be used to reahze a quantum gyroscope characterized by two important superfiuid 
effects: the reduced value of the inertia of the sample and the quantization of the angular 
momentum associated with the vortex. The proposed gyroscope is characterized by the 
precession of the symmetry axis of the condensate around the symmetry axis of the con- 
fining trap, the precession frequency being fixed by the angular momentum carried by the 
vortex line. The experimental realization of the proposed gyroscope would provide a further 
tool to explore the intriguing features exhibited by rotating Bose-Einstein condensed gases, 
including the stability [p!9| , ^ and the life time of vortex lines. 
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FIGURE CAPTION: 

Schematic picture of the BEG gyroscope. The condensate is suddenly rotated with respect 
to the symmetry axis of the trap (z-th axis) and exhibits precession in the presence of a 
quantized vortex line. 
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